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A new cluster model solution to the long standing nuclear structure problem of describing the 
anomalously long lifetime of 14 C is presented. Related beta-decay data for 14 to states in 14 N, 
gamma decay data between low lying positive parity states in 14 N and the elastic and inelastic 
magnetic dipole electron scattering from 14 N data are all shown to be very accurately described 
by the model. The shapes of the beta spectra for the A=14 system are also well reproduced 
by the model. The model invokes four-nucleon tetrahedral symmetric spatial correlations arising 
from three- and four-nucleon interactions which yields a high degree of SU (4) singlet structure for 
the clusters and a tetrahedral intrinsic shape for the doubly magic 16 O ground state. The large 
quadrupole moment of the 14 N ground state is obtained here for the first time and arises because 

t-H , of the almost 100% d-wave deuteron-like-hole cluster structure inherent in the model. 
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CN \ I. INTRODUCTION 

c : 

The history surrounding the existing explanations of the 14 C lifetime has involved many attempts (ll-fioj. but none 
of them have been completely satisfactory. As noted in the most recent publications [10| a full understanding of all the 
. data is expected to require three- nucleon interactions and/or clustering considerations. Unfortunately, until now, the 
details of clustering considerations and the nature of the three- (or more) nucleon interactions have not been addressed. 

i— i , Much of the previous work has involved A = 14 wave functions based on two p-shell holes in the closed shell reference 

i-^h ■ state (which is the 16 O ground state with {ls} 4 {lp} 12 structure). This model is still the one used in the most recent 
structure publication [ic|. Deviations from the simple two- hole shell model have Q invoked additional multi-particle 
/multi-hole states but the calculations show poor convergence and do not provide any accurate description of the 
beta decay data. Other investigators (Til . fl2l | chose to use phenomenological admixtures of the possible two-hole 
angular momentum configurations e.g., in L-S coupling: { 3 Si, 1 Pi, 3 Di} for the 14 N J p = 1 + , T = ground state and 
{^o^Po} for the J p = 0+,T = 1 isospin states in 14 C, 14 N and 14 0. These type of attempts were criticized (IHH 
I , for being inconsistent with the conventional strong j-j coupling shell model. As we shall show the phenomenological 
^ ■ approach appears to be closer to the model used here and with more searching might have resulted in the wave 

CO function admixtures presented here. 

In view of the failure of the shell model approach to provide a complete description of the beta decay and allied 
data in the A = 14 nuclei it does appear that the shell model approach is not the optimal choice as the basic picture 
for the 14 C beta decay problem. That this appears to be the reality of the situation led this investigator to invoke 
a more realistic model of the closed shell nucleus 16 which includes multi-nucleon correlations. Such a correlated 
model has a l ong er history than the shell model as it dates back to Wheeler [l4| in 1937. This was followed by other 
investigators [l5lfla| and for 16 O the alpha-particle cluster model relying on the similarity with the methane molecule 
CH4 was used to describe the energy levels as rotational vibrational states of a tetrahedral molecule in which the 
H-atoms were replaced by alpha-particles and the C-atom at the center was eliminated. All of the early work and 

Iv> j almost all of the later efforts with the alpha-particle model (see reference (T3| for a review) have assumed the alpha 
clusters to be uncorrelated {Is} 4 configurations as in the simple spherical shell model. The lack of correlations within 
each cluster and the assumption of spherical intrinsic states leads to difficulties in obtaining accurate predictions with 
the spherical alpha particle cluster model. In particular we know of no attempts to describe the lifetimes of 14 C or 
14 using such models. 

Quite early in the history of clustering this investigator proposed a model [18] which synthesized the simple cluster 
model with the shell model by introducing quark degrees of freedom into the bound states of nuclei with A = 2, 3 and 
4. The initial work (summarized in [19]) emphasized that many-body forces were to be expected and that the three 
and four nucleon bound states would have spatial correlations corresponding to point group symmetries D^h and 
for A — 3 and 4 respectively. A specific model [20j for the two-nucleon systems based on quark dynamics and one 
pion exchange gave a realistic description of the deuteron and the phase-shifts for low partial waves. The parameters 
of the model are consistent with the one nucleon non-relativistic quark model and the basic symmetries of QCD and 
chiral symmetry are adhered to. This approach was then extended consistently (2l1 ) to the spin independent part of 
the three-nucleon interaction which showed that an equilateral triangle configuration with D^ t point group symmetry 
for the nucleons was strongly favored. The strong repulsive interactions of up to 2 GeV between each pair of nucleons 
leads to a hole in the charge density distribution at the nuclear center as originally suggested by the authors of 
the experimental work [22T ]. 

The 4 He ground state is expected to have tetrahedral intrinsic spatial symmetry since it maximizes the three- 
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nucleon triangular configurations occurring on the four equivalent faces of the tetrahedron. Indeed as indicated in 
the elastic electron scattering data for 3 He and 4 He are very well described with equilibrium radii corresponding 
to triangular and tetrahedral geometry. Again the hole at the center of the 4 He charge distribution is well described 
by the model. Such spatial correlations are totally symmetric representations of the orbital angular momentum 
rotation group 0(3) provided that the intrinsic configuration is rotated through the three Euler angles with equal 
weight and no parity change under inversion. These spatial point group symmetries lead automatically to totally 
antisymmetric SU(A) states for spin S and isospin T, i.e., S = 1/2 = T for A = 3, and S = = T for the alpha 
particle. The 3 He, 3 H ground states are then simple one nucleon-hole states in the alpha-particle. Here we focus 
on the alpha particle which has J p = + and consequently L = only for its total orbital angular momentum (as 
also occurs in the simple {Is} 4 shell model configuration). At this point we note that the quark model discussed 
above showed significant quenching of the one-pion exchange tensor interaction as the nucleon- nucleon separation 
distance decreased and for heavier meson exchanges between quarks there was essentially no interaction. The 
spatial symmetry of the alpha- particle intrinsic state in this model is therefore expected to lead to very weak spin 
dependent contributions to the ground state cluster wave function. In what follows we will assume only the leading 
L = = S = T state to be present in the ground state of 4 He and in the intrinsic states of embedded four nucleon 
alpha-like clusters in 16 O. 

II. BASIS WAVE FUNCTIONS 

The cluster model wave functions used here show a strong resemblance to the shell model basis states and indeed 
for A — 14 the two- hole states have the same total orbital (L), spin (S) and total angular momentum J as those 
used in shell model states. For comparison purposes we use the same notation as that of Genz et al [ll| and in L-S 
coupling the most general wave functions are 

| 14 N, J p = 1+ T = 0) = a 3 S*i + 1 P 1 + 7 3 L>i (1) 
| 14 N*,J p =0+,T=l)=eo 1 5o+7 7o 3 P (2) 



| 14 C,J p =0 + ,T = l)=e- 1 5 +7 ? _ 3 F (3) 

| 14 O,J p =0+,T = l)=£+ 1 S + V+ 3 Po (4) 

wherein only the L p = 1 + (denoted by P) states have the same angular momentum substructure as the shell model 
uncorrelated p~ 2 hole states. In what follows the angular momentum substructure of the S and D cluster model states 
are not the same as the shell model uncorrelated p~ 2 hole states. The normalization of the above four configurations 
however is the same, i.e., 

a 2 +p 2 +1 2 =g+ V 2 = l (5) 

with l = 0,+and — . The values of £, b ,rj L allow for a possible isospin triplet symmetry breaking which is needed to 
describe the difference in the log(/.Ai) values for the /3 + and (3~ decays of 14 and 14 C leading to the ground state of 
14 N respectively. 

The reference state for the cluster model hole states is also the J p = + ,T = ground state of 16 O, but here 
we assume the reference state is highly correlated with four alpha-like intrinsic clusters with their centers of mass 
having equilibrium points at the four corners of a tetrahedron as in Fig. 9 of [19]. As in 4 He the ground state of 
16 O is found by rotating over all Euler angles with equal weight and requiring no change of parity under inversion of 
all the spatial coordinates. The strong cluster correlations between nucleons in the same alpha-like intrinsic cluster 
and the much weaker correlations between two nucleons in different intrinsic clusters leads to the assignment of 
S- and Z?-states for two nucleons taken from the same cluster and P-states when the two nucleons are taken from 
different clusters. In the cluster model presented here the individual clusters have SU(4) singlet structure and for 
four clusters satisfying identical boson symmetry then the 16 reference state will be pure L-S coupled. Removing 
a pair from an individual cluster allows the pair to have antisymmetric SU(4) quantum numbers, which yields only 
(S = 0,T = 1) and (S = 1,T = 0) SU(4) hole states. This in turn requires the relative motion of the two nucleons 
to be in an even parity spherical harmonic state. As discussed earlier the ^-value of the relative motion within 
a tetrahedral cluster is taken here to be zero. The total orbital angular momentum can be L — or 2 if the non- 
zero contribution comes from the motion of the center of mass of the two nucleons relative to the core. However 
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when the two nucleons are taken from different clusters the allowed SU (4) pair states must be symmetric in order to 
satisfy the bose permutation symmetry between the two identical clusters. In this case the allowed SU(4) states are 
(S = 0, T = 0) and (S = 1, T = 1) which requires their relative motion orbital angular momentum to be odd valued. 
The ground states of the mirror nuclei 15 N and 15 have J p = l/2~ and as in the shell model they are represented 
by proton and neutron p^-hole states in the ground state of 16 respectively. In the shell model the two p-holes can 
only have antisymmetric spatial states if L — 1 and this also holds true for the cluster model. Consequently we assign 
the P configurations to cluster model configurations in which individual p-holes are taken in different clusters. For 
the T = 1 states in A = 14 there are only two basis states: (a) the state involving a pure L = dinucleon cluster 
extracted from the same alpha-like cluster and (b) the 3 Po state involving a pure L = 1 with each p-hole taken from 
a different alpha-like cluster. Similarly for the T = basis states in 14 N there are two types of basis states: (a) the 
3 Si and 3 Di states are pure L = and 2 deuteron- like clusters extracted from the same alpha-like cluster and (b) 
the 1 P\ state with L = 1 only and with each p-hole coming from a different alpha cluster. 

The foregoing is important because it matters when considering the energy matrix for the individual systems. It 
also matters when considering observables that are dependent on the angular momentum substructure of the orbital 
L-states or on their radial wave functions. This is vital in the case of the quadrupole moment of the ground state of 
14 N and in describing elastic and inelastic electron scattering data. The /3-decay data, Ml gamma transitions in 14 N 
and the magnetic moment of the ground state of 1 N are independent of the substructures discussed above as these 
observables depend only on the admixture amplitudes defined in eqs. (Q][3]). 

III. ENERGY MATRICES 

For T = 1 states with two basis states the 2x2 matrices involve the coupling between the 1 So and 3 Po states which 
require a spin-orbit interaction which is antisymmetric in spin space and also in orbital space. Such an interaction can 
be constructed from the sum of one-body operators for each nucleoli and is a superposition of the nuclear spin-orbit 
VN so and the electromagnetic spin-orbit VE so , see j23|- Specifically these interactions for each nucleoli are 

VN so (neutron) = cr N ■ grad{UN(r N )} x p N = y7V so (proton) (6) 
VE so (neutron) = [a N ■ grad{UE(r N )} x p N ]fi n (7) 



VE so (pmton) = [* N ■ grad{UE{r N )} x p N }{n P - 1/2). (8) 

In these generalized spin-orbit interactions we assume that UN(r) and UE(r) are the nuclear and Coulomb potentials 
which can involve tetrahedral harmonics with orbital- values of 0, 3, 4, 6 etc. For convenience the nuclear terms in eq. 
© are taken to be the same for neutrons and protons and any differences (which should exist ) are taken to be included 
in the overall magnitude of the electromagnetic spin-orbit terms. The matrix elements for the T = 1 2 x 2 matrix we 
denote by Hg S , Hp P and Vg P (= V PS ), wherein the S and P subscripts imply the 1 ^ and 3 Po basis states. Values for 
these matrix elements are different for 14 C, 14 N (T — 1) and 14 and for diagonal elements one has the "unperturbed" 
energies of the system whereas the off-diagonal elements are simply the matrix elements of the spin-orbit interactions 
given above. The Vg P matrix elements are charge dependent because the nucleon magnetic moments /x ra and n P are 
of opposite sign. We characterize Vg P for each member of the isospin triplet for the nuclear spin-orbit matrix element 
by Wnuc and the strength for the electromagnetic spin-orbit for two proton holes by v c \. Consequently the values of the 
matrix elements Vg P for the isospin triplet are given by v nuc + f(N, N)v e \ with f(p,p) = 1, f(n, n) = fx n / (/x p — 1/2), 
f(p,n) = (f(p,p) + /(ft)ft))/2- The amplitudes ^ t ,ry t for the three values of t = — ,0 and +, corresponding to 14 C, 
14 N* and 14 respectively, are found by diagonalizing the three 2x2 matrices when specific values of the unperturbed 
energy differences Hp P — Hg S , v nuc and v e \ are used to fit the experimental data. In particular the experimental 
values of the three energy differences between the ground states of the isospin triplet and the corresponding first 
excited states with J p = + ,T = 1 are constraints to be satisfied by the eigenvalues of the three diagonalizations. 
These energy spacings are 6.589 MeV in 14 C, 6.305 MeV in 14 N* and 5.91 MeV in 14 0. This leaves two parameters 
out of the five input parameters for the T = 1 sector to be determined. 

For the J p = 1 + , T = case in 14 N there are three basis states which can be labeled by D, S and P corresponding 
to 3 Di, 3 Si, and 1 P\. We expect the Z?-state to be the least bound state in the cluster model as it is even in all shell 
models used historically In the cluster model we expect after diagonaliztion that the D-sX&te occupation will be close 
to 100% and will be the ground state of 14 N. The other two eigenstates should be the 1 + states at 3.948 MeV and 
6.204 MeV respectively. The diagonal matrix elements of the 3x3 energy matrix will have two unperturbed energy 
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input values e.g., Hg S — Hp D and H pp — H^ D and after diagonalization the eigenvalue of the system should have 
the energy splittings of 3.948 MeV and 6.204 MeV respectively. The other three input values are the matrix elements 
Vjj P , Vg P (which are non-zero from the spin-orbit interaction in a similar manner to the T = 1, Vg P discussed above) 
and Vjj S with the latter expected to be very weak since it can only arise from a scalar product of rank two tensors in 
spin and orbital spaces. This expectation is based on the fact that tetrahedral symmetry does not have quadrupole 
harmonics since this suppresses all tensor interactions in the 16 system. The T = sector (like the T = 1 sector) 
involves five input parameters and fits to the energy differences reduces this number of parameters to three which 
must be determined from other data. Thus overall we have five variables in the energetics, two from the T = 1 sector 
and three from the T = sector, and at first sight these could be determined by the three /3-decay data sets for 
14 C(/3~) and 14 0(/3 + ) going to the ground and first excited J p = 1 + states in 14 N. The five observables are the three 
\og(fAt) values and the shapes of the j3~ and f3 + decays to the 14 N ground state. The shape of the (3 + decay to the 
3.948 MeV state in 14 N is not measured and in any event is expected to be constant as this transition has a log(/^t) 
of 3.138 corresponding to an unhindered Gamow- Teller transition. Unfortunately, as shown by Towner and Hardy 
[24j . one needs to invoke renormalized axial (gA) and magnetic (gi and g s ) couplings to obtain the correct results for 
the strongly hindered j3- decay data and the radiative width of the T = 1 N (2.313 MeV) state. 

IV. RENORMALIZED OPERATORS 

The operators needed in the remainder of this paper are the free nucleon coupling constants gA,gip,gimgsp and 
g sn corresponding to the axial vector (gA = 1.2695), the orbital g-factors (gi p = 1 and gi n — 0) and spin g- factors 
(g sp /2 = 2.79285 and g sn /2 = —1.91304), with all the magnetic couplings being in units of nuclear magnetons (n.m.). 
The renormalized gA for the j3- decay studies in A — 14 is taken from the 15 0(/? + ) 15 N mirror state transition as 
suggested by Towner and Hardy. The transition is assumed to take place between single p-shell holes in the 16 O 
reference state and results in a renormalized value of g* A — 1.0885 = (gA — 0.181) which yields the log(f^t) = 3.644 
for the Gamow- Teller component [25[ of the 15 0(/3 + ) 15 N transition. The renormalized magnetic operators are taken 
here initially to fit the magnetic moments of the ground states of the mirror nuclei 15 N and 15 based on a single 
p-orbital-hole in the 16 O ground state. In looking at these magnetic moments using the simple formulas for a p-hole 
with j p = l/2~ given by 

Hi = l/3(2g H _g si /2) = g Jt /2 (9) 

where for i =free neutron g n i = and for a free proton g p i — 1 and similarly g ns /2 = —1.91304 n.m. and g ps /2 = 
2.79285 n.m. one finds values of the moments g^/2 as -0.26428 n.m. for 15 N and +0.63768 n.m. for 15 0. These 
are not in good agreement with the experimental values of -0.28319 n.m. for 15 N and +0.7189 n.m. for 15 0. It is 
necessary to use renormalized magnetic couplings as pointed out by Towner and Hardy [24j . They obtained these by 
including bound state shell model corrections arising from core polarization and meson-exchange currents. 

In the cluster approach the p-hole arises by breaking the individual alpha-like clusters in which it is embedded. 
This suggests that an initial guess for the renormalized 3^/2 values in the A=15 systems should be given by 

g* t (A = 15)/2 = g 3l (A = 15)/2 x {g* t (A = 3)/ gji (A = 3)} (10) 

in which g^(A = 3)/2 are the observed moments for the s-hole states in a free alpha particle which are -2.12750 n.m. 
and +2.97896 n.m. for the neutron and proton s-holes respectively. The values of gji(A = 3)/2 are taken to be the 
free nucleon magnetic moments. The renormalized magnetic moments fi* for A = 15 are found to be -0.28189 n.m. 
and +0.7092 n.m. for 15 N and 15 respectively. This cluster model approach to the renormalization of the magnetic 
moments of the nucleon- holes is remarkably accurate in obtaining values of the A = 15 magnetic moments which are 
within 0.5% and 1.4% of the observed values for 15 N and 15 respectively. 

The above discussion for finding renormalized magnetic coupling constants that describe the A — 15 mirror states 
appears to be another validation of the alpha-like cluster model for the 16 reference state. However further mod- 
ifications to the values of <?*, and g* si are needed to obtain a consistent and completely accurate description of the 
magnetic moment data for the ground state of 14 N as well as the A=15 mirror pair states. In using the result from 
(ITU1) above we infer that gf is approximately 1.1 , g* n is approximately 0.0, \i* sp is approximately 3.0 and n* sn is close 
to -2.15 (or if gl n = then = —2.1567 so that the magnetic moment of 15 is exactly reproduced). Best results 
are achieved with the values: = 1.112, g\ n = 0, [i* sp = +3.0735664 and ^ sn — —2.1567 using the wavefunctions 
that fit the (3 - decay data as discussed below. 

Not only does this set of renormalized couplings fit the magnetic moments of 14 N (/i — +0.403761 n.m.), 15 N and 
15 but also exactly fits the magnetic moment of the first J p — 3~ state in 16 0. The measured value for this 3~ 
state is fi = +1.668 n.m. corresponding to gj with the isoscalar gyromagnetic g = +0.556 (error is .005) . In the 
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tetrahedral model this 3~ state is a collective rotational excitation of the + ground state of 16 for which the g 
factor is Z/A — 1/2 if the orbital gu is 1 for protons and for neutrons. Using the cluster renormalized orbital g- 
factors given above yields g = g* p /2 = +0.556 in perfect agreement with the data. In obtaining the magnetic moment 
for the ground state of 14 N the shell model formula [llj 

H = 2- 1/2 {(/i P + li n )W 1 /2- 1 / 2 + W 3 } (11) 

(in which W\, W3 are (2a 2 — j 2 ), ((3 2 + 3j 2 /2)/2~ 1 / 2 respectively), is modified for the cluster model to fx* given by 

^ = 2-i/ 2 {^ d W 1 /2- 1 / 2 + g: p W 3 } (12) 

with n* d being the renormalized isoscalar magnetic moment of the deuteron-like hole state in 14 N. The deuteron-hole 
can involve internal orbital angular momentum of / = and 2 due to the two-body tensor interaction between the 
neutron- and proton-holes in the same alpha-like cluster. We use the equivalent parameterizations for = /d(/-t p +/-tn) 
or n* d = f d ((J.p + (J-n) the only new parameter is fd (or equivalently f d ) since we use the renormalized nucleon-hole 
magnetic moments given above. If f d is unity then the probability Pd of any I = 2 state in the deuteron-hole is zero. 
For a small value of Pd we should have /J being slightly less than unity. Using the renormalized operators which 
fit the A = 15 and A = 16 magnetic moments then the 14 N moment given by (fT2"j) is 0.382 n.m. if /? = 1 and is 
0.403761 n.m. when we choose /J = .950706 (or fd — .990755) which is consistent with a small value for Pd- The 
value of Pd — 8.3% obtained here using 

P D = (2/3)(l - /,*)(,;; + + (x* n - g* lp /2) (13) 

is not inconsistent with any information concerning D-states in the alpha particle. However by allowing for 1% errors 
in the magnetic moments of 15 and 16 one can obtain fits to the magnetic moment of 14 N which have values of Pd 
as low as 4.7%. 



V. RESULTS FOR BETA-DECAY OBSERVABLES 



For Gamow- Teller (GT) transitions in the A = 14 nuclei we restrict our considerations to J p = + initial states 
leading to J p = 1 + states in 14 N. The GT matrix element is given by [24] 

MGT^e 1 / 2 ^-^ 172 ) (14) 

where we include the renormalized axial coupling g* A in the definition of MGT. With this definition the f^t is given 
by 

f A t = 6146/|MGT| 2 s. (15) 

in which the Ja rate functions are corrected to include the effects of the nuclear structures via the shape function 
C(Z, W) used by Towner and Hardy 

C(Z, W) = |MGT| 2 fc(l + aW + fiijib/W + cW 2 ) (16) 

following the format of Genz et al that involves the a, b, c variables which are dependent on the details of the nuclear 
structure model. To the accuracy needed for A — 14 beta-decay shapes we use /ii7i = 1. The parameter k is used 
to fit the data and W is the total electron (positron) energy in units of the electron rest mass energy. The data [27| 
for the 14 C(/3 - ) 14 N decay gives information on the slope parameter (a) whereas the data [28[ for the 14 0(/3 + ) 14 N 
decay gives information on all three parameters (a, b and c) . The expressions for a, b and c as well as the nuclear 
matrix elements are given in detail in Genz et al [Tl| and also in Garcia and Brown (2(| who noted that the term 
denoted by V4 has the opposite sign from that given by Genz et al. Our calculations use the sign choice of Garcia 
and Brown for V4 as it is consistent with the relation between 2-hole states and 2-particle states. We choose not to 
include these lengthy relationships as they are readily available in [Til [26| . It is important to note that the shape 
functions calculated using the formulas in Genz et al have some terms which are model dependent. In particular 
we have used renormalized g* for axial and magnetic couplings in place of the free nucleon g's and also replaced the 
oscillator length (b = 1.7fm in Genz et al) by the cluster model b c value obtained from the cluster model fit (using 
linear combinations of p-state radial harmonic oscillator states) to the inelastic electron scattering data as discussed 
in VI below. The shape function is sensitive to the choice of b c . The behavior of the slope parameters a_ and a+ for 
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Figure 1: The slope parameters a_ and a+ used in (|16[) for C(Z, W) for the 14 C(/3~) and 14 0(/3 + ) decays to the ground state 
of 14 N respectively as a function of the average oscillator length b c . The log(/At) values are independent of b c and correspond 
to the tabulated wavefunction admixtures given in section V. 

14 C(/3~) and 14 0(/3 + ) decaying to the 14 N ground state using the wavefunctions shown below are shown in Fig. [T]as 
a function of the average oscillator parameter which is denoted by b c in Fermi units. 

The cluster model wavefunctions that fit all the data in terms of the coefficients a, /?, 7, £ t , rj L are given for 
T = 0, J p = 1+and for T = 1, J p = 1+ in Table Hand Table HI respectively. 

a ft 7 

| 14 N,£ = 0) .0169003 .1860092 .9824026 

14 N,£ = 3.948) .7600690 -.6407632 .1082473 

14 N,£ = 6.204) .6496225 +.7448645 -.1522089 

Table I: Numerical values of a, /3, 7 for the three 1 + states in 14 N which result from the diagonalization of the symmetric 3x3 
matrix as discussed in the text. 



\ 14 C,E 


= 0) 


.9865400 


.1635201 


| 14 N,£ = 


2.313) 


.9909549 


.1341954 


\ 14 0,E 


= 0) 


.9945400 


.1043556 


\ 14 C,E = 


6.589) 


-.1635201 


.9865400 


| 14 N,£ = 


8.616) 


-.1341954 


.9909549 


| 14 0,B = 


5.910) 


-.1043556 


.9945400 



Table II: Numerical values of £ t , r/ L for the two + states in each of the isospin triplet of A — 14 nuclei which result from the 
diagonalization of the symmetric 2x2 matrix as discussed in the text. 



These structure coefficients are the eigenstates obtained by diagonalizing the symmetric 3x3 matrix for T — with 
off-diagonal matrix elements (in MeV) :Vg D = -0.28864, Vg D = -0.9772 , V<$ P = 1.079353 and diagonal unperturbed 
energies Ess = 4.7088, Epp = 4.873 and Edd = 0. Similarly there are three symmetric 2x2 matrices for the T = 1 
systems and their matrices have elements given by Table Hill in which the unperturbed energies have an average value 
for the isospin triplet which is 6.03 MeV, that is remarkably similar to the energy splitting of the first two J p = + 
states in the reference system of 16 0. 
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Vsp Ess Epp 
C -1.062993 6.237 
N -.8381848 6.076 
O -.6133766 5.7813 

Table III: Numerical values of the elements of the 2x2 energy matrix for the isospin triplet of A = 14 nuclei. 



The major effect of the electromagnetic spin-orbit in the T = 1 isospin triplet states lies in the relatively large 
differences of the r\ L coefficients. In particular the magnitude of 77+ for the 14 ground state is about 36% smaller 
than rj_ for the mirror state in 14 C. Indeed it is this large difference in the cluster model r\ coefficients which leads to 
an explanation of the large ratio of the MGT elements for 14 0(/3 + )/ 14 C(/3~). It is vital to understand that it is the 
cluster correlations between the nucleons that damps the nuclear spin-orbit (v auc ) and enhances the electromagnetic 
spin-orbit (v c \). The values used here are v nuc = —0.817874 MeV and v c \ = —0.245119 MeV. As defined above 
u e i is appropriate for the two proton-hole 14 C system and the values of Vg P are given using the superpositions of 
Unuc + f(N, N)v e i as given in III above. Although the symmetry breaking in 77+ and 77_ is large the overlap of the 
14 C and 14 ground states is 0.99822 which involves an overall symmetry breaking of less than 0.2%. 

Using the above wave functions and b c = 2.25 fm leads to the values of MGT, log(/^i) and the shape parameters 
a,b,c given in Table IIVI 



Model 



MGT 



log(f A t) k 



14 0(£+) 14 N(£ = 0) 
Cluster 
Cluster 
PBWT 
Expt: 14 Q 



0.01493 7.440 

0.01373 7.513 

0.01480 7.448 

0.018(0) 7.284(7) 



1.600 -0.056 0.033 0.002 
1.990 -0.060 0.036 0.003 
1.605 -0.054 0.035 0.002 



14 0(/?+) 14 N(£; = 3 
Cluster 

.14 



Expt: i4 



2.119 3.136 1.0 
2.119(39) 3.138(16) 



0.0 0.0 0.0 



14 C(/T) 14 N(£ = 0) 
Cluster 
Cluster 
PBWT 
Expt: 14 C 



-0.00237 9.040 

-0.00343 8.718 

-0.00343 8.718 

0.002(0) 9.040(3) 



-0.231 0.125 0.023 
-0.176 0.063 -0.003 
0.607 -0.235 0.013 0.013 

-0.23(2) (100-160 keV) 
-0.17 (50-160 keV) 



Table IV: P + and /3~ transition results for MGT, log(/ytt) and k, a, b, c shape parameter values using cluster wave functions 
from the text and PWBT shell model results (with their renormalized operators) from Towner and Hardy (24| . Expt. values 
for |MGT| and log(fAt) are from |25j |. The slope parameter a for 14 C is taken from Table II of [27l | using the fitted value of 
a = —0.45(4) MeV -1 in the energy range 100 — 160 keV. In electron rest mass units this becomes a = —0.23(2). For the wider 
energy range a in electron rest mass units is —0.17. 



For the special case of 14 0(/3 + ) 14 N(_E = 0) the present work yields almost identical results to those of Towner and 
Hardy (who are the only ones to use renormalized couplings) . The fit to the experimental data they obtained for the 
shape function using the PWBT [2j| shell model basis shown in Fig. [5] has a slightly larger % 2 than that obtained 
with the cluster model. The best fit to the data is obtained using a small change in r] + to 0.1085 which gives the 
dashed curve in Fig. [5] and corresponds to a log(/^t) of 7.51 and an MGT of 0.01373 for the 14 0(/3 + ) transition. The 
sensitivity of the slopes a + and a_ and the corresponding log( f At) values to the values of 77+ and 77_ is shown in Fig. 
[3J The sensitivity of a+ to r/ + is very much less than the sensitivity of a_ to r\- which in the cluster model is very 
sensitive. 

The cluster 14 C(/3~) 14 N(_E = 0) results are significantly different from those of Towner and Hardy as our wave- 
functions give the usual value of 9.04 for the log(/^t) when the slope parameter a_ = —0.231 rather than the PWBT 
shell model log(/>it) of 8.72. We have been unable to exactly pin down the source of this difference. We note that the 
main thrust of the Towner and Hardy [24| investigation was to show that CVC could be satisfied accurately by using 
renormalized g factors and in the case of 14 0(/3 + ) 14 N gs by adjusting the PWBT ground state of 14 N they obtained a 
good description of the /3 + shape function. We are concerned however that the Towner and Hardy calculations used 
a different wave] 'unction for 1A N{E = 0) in the 14 C decay from the one used in the 14 decay and recently Negret 
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1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 3.9 

Positron Kinetic Energy (MeV) 

Figure 2: The shape function C(Z, W) for the 14 0(/3+}i 4 N(£ = 0) transition as a function of the positron kinetic energy 
(MeV). The data points with error bars are taken from [2|| including corrections given in [24l |. 
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11- T1 + 

Figure 3: The values of a_ and a+ as a function of the respective initial state admixture coefficients denoted by tj_ and 77+ 
are shown by the continuous lines in the left and right side panels respectively. The corresponding values of the log(/Ai) are 
shown by the dashed lines and for the 14 C(/3~) case shows the strong variation in a_ and log(/ J 4*) for small changes in r\- near 
0.1635201. 



et al [30] implied that Towner and Hardy had been able to account for the large asymmetry in the mirror log(ft) 
values. Physically these mirror decays have the same final state in 14 iV and one cannot explain the large ratio of their 
MGT values by simply modifying the final state. Indeed Towner and Hardy did not claim |24| to have explained the 
large asymmetry in these mirror \og(ft) values; their focus was on reconciling the shape- correction form factor 
with the Ml matrix element in 14 7V. As pointed out by early workers 0, H[ the symmetry breaking of these mirror 
transitions must arise from symmetry breaking interactions in the initial mirror states. Here we have specified that it 
is primarily the electromagnetic spin-orbit interaction in the initial states that causes the large asymmetry. It is also 
worth noting the opposite signs of the MGT values for these mirror transitions because the shape parameters do not 
come out correctly unless the mirror MGT values have opposite signs. This is most readily seen in Fig. [3] where the 
values of ?7_ between 0.14 and 0.15 can yield log(/yii) values between 8.5 and 9.4 with MGT > but also yield a 
positive slope for a_ in complete contradiction to experimental data. 

One of the difficulties remaining is the uncertainties in the experimental data [27] for the 14 C(/3~) transition. In 
particular the value of the slope parameter depends strongly on the range of electron energies used as indicated in 
Table IIVI Also the accuracy of the data apparently did not allow any information to be determined for the b or c 
coefficients in C(Z,W) so that a linear form Cl(Z,W) = (1 + a*W) was used [27J to extract the effective slope a* 
for each energy range. The value of a* = —0.45(4) MeV -1 for the 100-160 keV range is apparently believed to be the 
favored value for a* in ref. [13. In the case of the shell model calculations given in ref. |24j one can least squares 
fit the models shape functions over the energy range 100-160 keV with the linear form to obtain the effective a*. We 
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find a* = -0.203, -0.214, and -0.195 for the three models labeled CK, PWBT, and MK respectively in Table III 
of Towner and Hardy. The corresponding values of a are —0.215,— .235 and —.207 which are only 6% — 9% higher 
than their respective a* values. For the cluster model with log( f A t) = 9.040 and a — —0.231 a linear fit to C(Z, W) 
gives a* = —0.232 because of the strong cancellation of the b and c terms in this case. The cluster case where the 
log(f A t) = 8.718 and a = —0.176 when linearized yields a* — —0.207 which is consistent with the shell model values 
for a*. If one accepts the best value for a* = —0.23 then the cluster wavefunction with log(/^t) = 9.040 is favored 
over all the other models. Apparently only more accurate data for the 14 C(/3~) transition can provide the necessary 
information on a, b and c, or a*. 

The /3+ decay of 14 leadingto the 14 N (E = 3.948) state calculated here has a log( f A t) of 3.14 and an MGT 
of 2.119 in perfect agreement |25| with experiment. This is not discussed by anyone else except for Genz et al (TlT | 
who quoted that their model yielded a log(f A t) value of 2.87. In our case this latter transition was part of the fitting 
procedure whereas Genz et al did not include this transition as part of their fitting procedure. However we believe 
that calculations of this transition to the first excited 1 + state in 14 N are an additional restraint that all models should 
include. 

The major reason that renormalized couplings were used is to understand whether the radiative (-Ml) width of the 
14 N,T = 1 state at E = 2.313 is consistent with the model wave functions. Again the model states above yield a 
value of T 7 = 6.7meV but only if renormalized values of the isovector coupling constants from above are used thereby 
satisfying the conserved vector current requirements. The formulas are given by Garcia and Brown [26] in their eq.(25) 
in terms of the structure coefficients V\ and V3 given in [Tl1 | and with our renormalized magnetic couplings one needs 
in order to obtain T 7 that 

l 2 ~ 1/2 (/^ - Mn)^i(7) + 9i p v s\ = 0.256 n.m. (17) 

with an associated error of 0.006 n.m. The value obtained here is 0.256 n.m. since the radiative width T 7 of 6.7(3) 
meV was part of our fitting procedure for the coefficients a, ft, 7, £0 and 770. The recent work of Holt et al [T3| also 
considered the recent experiments [30( that determined BGT = |MGT| 2 /(2 Ji + 1) values from the 14 N ground state to 
excited states of 14 C and 14 using the charge exchange reactions 14 N(d 2 , He) 14 C and 14 N( 3 He,t) 14 0. For the three 
final states in each of 14 C and 14 labeled as 0^~,0j we have wavefunctions as listed above and for th e it states at 
11.31 MeV and 11.24 MeV respectively we use a pure 3 Pi two nucleon-hole in 16 as did Amos et al [13] . For the 
0^~ states the beta decay data is more reliable and we already fitted their BGT values. For the 0^ states we obtain 
BGT = 0.028 for both these transitions which are similar to the very small values shown in Negret et al [3(j. By 
using the definition in Holt et al for B(GT) corresponding to the inverse transition (note that g A is not included in 
their definition of B(GT)) we obtain for the 14 C case that B(GT) = .071 which appears to agree closely with the 
experimental value shown in their Fig. [5] for this transition. The Holt et al theory result for this transition is at 
least three times too large even with the modified tensor interaction included. The Ml transition from the ground 
state of 14 C to the 1+ state at 11.3 MeV in 14 C was first discussed in [l3j as a way to obtain information on the 3 Po 
component in the ground state of 14 C. The radiative width T 7 of this 1^" state was measured [3lJ to be 6.8 meV by 
extrapolating the inelastic electron scattering from 14 C data to the photon point However such extrapolations can 
be quite inaccurate and it is more reliable for this if state to consider the BGT measurements of Negret et al. The 
calculated value of BGT to this state in the cluster model is BGT = 0.082 or B(GT) = 0.069 which appears to be 
only about 30% below the experimental value. 

The situation for the transitions to J p = 2 + states from the 14 N ground state is beyond the scope of this work 
as it is too complicated for us to calculate with any accuracy in the cluster model since there are three 2 + states 
observed in each mirror system with significant strength and only two 2 + states with the simple two-hole structure in 
the 16 ground state. However the two simple states , 1 Z?2 and 3 P2, will have a strong BGT only for the 1 Di because 
the initial state has over 98% in amplitude in the 3 £>i. This means that the transitions to the three mixed 2 + states 
observed will dominate the GT transitions as suggested by Aroua et al Q and observed in the experiments of Negret 
et al. 



VI. ELASTIC AND INELASTIC Ml ELECTRON SCATTERING 

In this section it is necessary to go beyond the admixture coefficients a, /3, 7, £ t , r\ L and the angular momentum 
quantum numbers because the form factors for Ml transitions involve momentum transfer (q) dependence and four 
types of multipoles for each nucleon [32j. In the p~ 2 simple shell model it is customary to use a single p-shell harmonic 
oscillator wave function for each nucleon and this yields very simple spherical Bessel transforms for L = and L = 2 
amplitudes. These are given by Willey [32j as radial integrals over jj~(<ir) and a unit normalized lp radial density : 



0'o)ip,ip = (1 - 2x/3)e x , (j 2 )ip,i P = (2x/3)e 



(18) 
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in which x = q 2 b 2 /4 and b is the usual three dimensional harmonic oscillator length. In the case of 15 N there is elastic 
electron scattering data available [33| for the Ml form factor which we denote as -FY(g). In the simple shell model 
this Ml form factor involves a lp\/ 2 proton- hole description and the form factor is given by: 

F T {q) = Kaq{{! - 2x/3) + Ax/9{g pl + fi p )/ a} e - x F SN F c . m . (19) 

in which F$n is the nucleon size form factor and F c/m . for harmonic oscillator states is a simple Gaussian so that 
eT x F c m is replaced by eT x \ A ~ x )l A with A — 15. In this work we use a simple dipole form for F$n = (1 + 
.054675g 2 ) -2 corresponding to an rms radius of 0.81 fm. The observed magnetic moment fi is —0.2831888 n.m. and 
K = 2- 1 / 2 hc/{2nMc 2 Z) has the value 0.02124 fm for Z = 7. The above formula has been checked for the shell 
model calculations with various b values and reproduces the results given in [33] using free nucleon g p i, u p values. As 
Singhal et al [33| point out the simple shell model calculations overestimate the peak value of l-Frl 2 by 20-30% and 
underestimate it for q e ff beyond 2.4 fm -1 by large factors. Using our renormalized g*i, /i* in place of the free values 
in (|T9l) makes the shell model calculated \Ft\ 2 values even larger and makes no significant improvement in the large 
q behavior. Only by including configurations from the 2^/2 shell can the data be fitted |33| for all q values. 

In the cluster model for 15 N we expect the radial distribution of the proton-hole in the highly correlated tetrahedral 
alpha-like cluster 16 reference state to be considerably different from the lp independent particle radial wave function. 
Lacking a detailed theory for the many-body interactions it is convenient to use a simple expansion for the Fourier- 
Bessel densities (|18[) as a linear superposition of lp densities with different values of the oscillator parameter: 

F c . m .(j )i P ,i P = n- 1 ^(1 - 2x(n)/3)e-( A - 1 ^/ A , F c . m (j 2 ) lpAp = 2x{n)/3e^ A -^)/ A (20) 

n 

in which x(n) — q 2 b 2 j4: and n is limited to four. The four values of b n are found by fitting the Ml elastic electron 
scattering data for 15 N and 14 N. For the 14 N case we use the similar expressions to those used in [34] 

F elT (q) = qn- 1 ]T e-^- 1 ^' A {A + A lX (n))F SN (21) 



F mT {q) = q^ 1 £ e -M'WM( ft + B 1 x(n))F SN (22) 

n 

with the Bi amplitudes being given in terms of the structure amplitudes in L-S coupling as 

A = (2/3y/ 2 nK, A 1 = -(2/3) 1 / 2 A>*(PF 1 - W 2 )/3 (23) 

B Q = -(3/2) 1 /2 A -{2- 1 /2 (/i ; _ u* n )Vi + gt p V 3 }, B 1 = (2/3)V*K(£ - ^(V, - V 2 )/3 (24) 

in which /x is the magnetic moment of 14 N (calculated here to be the observed value of —0.2831888 n.m.). The 
coefficients Wi are defined in Genz et al [TTjj and in our model they are: 

Wi = 2a 2 - 7 2 , W 2 = -(4/5) 1/2 a 7 + (27/10) 1/2 /3 7 + 7 2 (25) 

in which the W 2 differs from theirs by the coefficient of 7 2 due to the deuteron-like hole in the same alpha-like 
cluster having the orbital angular momentum (L = 2) entirely in the coordinate connecting the center of mass of the 
hole pair to the center of mass of the reference system. Note that the quadrupole moment Q of 14 N is given here by: 
Q = (r 2 )/5{(16/5) 1 / 2 a 7 -/3 2 + 7 2 } which differs from the expression for Q in Genz et al only by our using unity instead 
of 7/10 as we also used in (fTT)) . The coefficients Vi are identical to those in Genz et al; V\ = — 2 1 / 2 (£oa: — ^o/?/3 1//2 ), 
V 2 = -((2/5) 1 / 2 Co 7 + 6^ 1/2 ?7o/3 + (9/20) 1 / 2 77o 7 ) and V 3 = - (2/3) 1/2 £o/?+ (2/9) 1/2 r l0 a - (5/18) 1 / 2 r ]0 7 which when used 
in pT|) yield with the wave functions in section V the result of 0.256 n.m. as needed for the radiative width to be 
6.7 meV. The value of the radiative width being 6.7 meV in the model used here is independent of up to 30% variations 
in the value of 770 . In this model the value of 770 must satisfy ?y_ > 770 > ?7+ and in general r]Q lies approximately 
half-way between r7_ and 77+ . Since the model is close to L-S coupling for all ground state wave functions and all the 
r/s are less than +0.2 to yield good descriptions of the beta-decay data to the ground state of 14 N it follows that the 
cluster picture with consistent renormalized operators g* is in full agreement with the requirements of CVC. 

The form factor for elastic electron Ml scattering has been fitted using the set of four b n values (all in fm units) 
bi = 2.85,62 = 1.95,63 = 1.82 and 64 = 1.32 and shows excellent agreement with the experimental data [34j in 
Fig. UJ The same set of oscillator lengths is used to calculate the elastic electron Ml scattering form factor for 
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Figure 4: The transverse magnetic elastic form factor for 15 N calculated using the cluster model parameters described in section 
VI and shown by the continuous curve is compared to the Bates data [33| . 

14 N and is compared with the experimental data in the lower panel of Fig. O The fit is of the same quality as 
that shown for 15 N in Fig. 0] The consistency of the model for A = 14 and A = 15 gives some credence to the 
use of the cluster modified lp density given in (20) above. The rms proton-hole radius rmsph is calculated using 
rmsph = {5{A - l)/8(6f + b\ + b\ + bl)/A + (0.81) 2 } 1 / 2 = 3.250 fm (A = 15), or 3.242 fm (A = 14) and it is the 
A = 14 value we need for the calculation of Q as defined above. Using the values of a, (3, and 7 given in section V 
above we find that Q = 20.18 mb that agrees well with the recent experimental values of 19.3(8)mb and 20.01(10) mb 
from the recent compilation (35|. A small correction to the calculated value of Q arises because of the small value of 
Pd needed to fit the magnetic moment of 14 N. The estimate of this correction is uncertain but one obtains a value of 
— Q^7 2 /10 in which Q* d is the quadrupole moment of the bound deuteron in 16 0. We anticipate that Q* d is most likely 
less than the Qd of the free deuteron (+2.86 mb) so that the correction is expected to lie between and —0.3 mb 
where the minus sign is because of the hole nature of the deuteron. Thus our estimate of Q for 14 N is 20.03(15) mb 
where the error is comparable to the error in the experimental value for Q. It is interesting to note that the early 
value of Q given in 1955 by Sherr et al was about 7 mb and it grew steadily over the next 38 years to the most 
accurate value of 20.01(10) mb in 1993 . 

It is also possible to use the simple relation used in [36j between the 16 charge radius and the 15 N charge radius 
for our model state for the ground state of 15 N: 

(r 2 )T ( 15N ) = { 8 /<V)ch( 16 0) - (rmsph) 2 /7} 1/2 = {8/7(2.71) 2 - (3.25) 2 /7} 1/2 = 2.62 fm (26) 

which agrees very well with the result given for the charge radius of 15 N in [36]. 

The inelastic electron scattering from 14 N leading to the first excited state at E = 2.313 MeV shows much more 
deviation from the shell model calculations than the elastic data does and because this transition involves a change 
of isospin from T = to T = 1 it has the connection to beta-decay as indicated in earlier work concerned with CVC. 
As we noted earlier for the beta-decay C(Z, W) shapes we need a knowledge of the average oscillator length squared 
which we denote by 6 2 and which we can assume should be the average of a set of 6 2 that fit the inelastic electron 
scattering form factor. Our fit to the data using (all in fm units) b\ = 3.512,^2 = 2.12,63 = 1-50, and 64 = 1.09 
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Figure 5: Measured elastic and inelastic Ml form factors for the ground state and 2.313 MeV transition in 14 N taken from [34 
(and references therein) are compared with cluster model calculations using the parameters defined in section VI. 



describes the data very well including the large q values. The value of b c from b\ = (b\ -\ 
which was used to calculate the shapes for the beta decays discussed in section V above. 



&D/4 is 2.25 fm 



VII. GAMMA DECAYS IN 14 N 



In most of the work on beta decay in the A — 14 system only the Ml decay of the first excited state has been 
calculated. In this section we focus on the gamma decay of the J p = 1 + at E = 3.948 MeV excitation. The gamma 
decay to the ground state involves two multipoles corresponding to E2 and Ml. The wave functions in section V 
gives the values of 0.0026 eV and 0.00042 eV for the radiative widths for E2 and Ml respectively which agree well 
with the corresponding experimental values 0.003 eV and 0.0004 eV from the TUNL compilation [37]. The gamma 
decay of the 3.948 MeV state to the T = 1, + state (2.313 MeV) is the almost completely dominant decay mode and 
is pure Ml. The wave functions in section V give a radiative width of 0.155 eV which is in good agreement with 
the experimental value of 0.091(30) eV. Overall the cluster model calculated gamma decay widths of the 3.948 MeV 
state are in satisfactory agreement with experiment but much of the gamma decay data in 14 N is quite old and new 
measurements could provide more accuracy on the transitions between the low-lying positive parity states. 



VIII. DISCUSSION AND CONCLUSIONS 



The results given in the preceding sections represent a convincing argument that the standard simple shell model is 
not an optimal starting point for describing the ground state of 16 O which, in turn, means it is not an optimal starting 
point for describing the low-lying states of A = 15 and A = 14 nuclei. In particular the cluster model need for the + 
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ground states of A = 14 nuclei with T — 1 to be more than 97% in the 1 S'o configuration (corresponding in j-j coupling 
to 1/3 probability for the p^J 2 component and 2/3 probability for the P7? 2 component ) contradicts the expectations of 
the strong spin-orbit j-j shell model ideology as pointed out by Talmi [8| ■ More important however is that the strong 
correlations leading to alpha-like clustering with T c i point group symmetry is that 16 is a tetrahedrally deformed 
nucleus and not spherical, as all shell model calculations use as a basic starting point. It is strongly deformed as the 
ground state rotational band [l!| with the sequence + ,3~,4 + ,6 + , has strongly enhanced BE3, BEA values for the 
transitions from the 3 — , 4 + states respectively to the ground state that are typical of a simple tetrahedral rotor model. 
There is no quadrupole deformation in the tetrahedral model of the ground state rotational band of 16 O as it violates 
the boson symmetry for four identical alpha-like clusters (tetrons?). We believe for N = Z even-even nuclei that the 
clustering dynamics of (N + Z) /4 identical tetron clusters determines the various multipole deformations with L > 1 
for each nucleus. 

We note that the next alpha-like cluster system to not have any quadrupole deformation in the lowest energy 
intrinsic state is 40 Ca as it appears to also be tetrahedrally deformed rather than spherical. In particular the lowest 
lying excited states are + {E = 3.35 MeV) and 3~(E = 3.74 MeV) with the BE3 value being 31 (W.U). Surprisingly 
at first is the fact that the binding energy of the last neutron (15.643 MeV) in 40 Ca is almost identical to the binding 
energy (15.664 MeV) of the last neutron in 16 . The difference is 21 keV and we do not for one moment believe this 
is an accidental degeneracy. We conjecture that the last neutron taken from one of the outermost tetron clusters will 
have three neighboring tetron clusters which are in the same close packed configuration as the four tetrons in 16 O. 
The remaining six tetrons in a tetrahedral 40 Ca are all spatially removed from the four containing the last neutron so 
that there is almost no interaction between the six spectators and the neutron being taken out. In short the similarity 
in the neutron binding energy in the tetrahedral model arises because the weakest bound neutron interacts only with 
the nearest neighbor clusters which is the same in 16 and 40 Ca. Of course this does not happen for the binding 
energy of the last proton because it sees the long range Coulomb interaction from all the spectators. 

In summary the new solution to the A = 14 beta-decay puzzle uses a highly correlated model in which the 
unperturbed LSJT basis has a simple orbital angular momentum selection rule which forbids the ground state GT 
transitions in the mirror /3-decays. This rule arises in this model because the L = 2, S = 1, J = 1, T = state is the 
most appropriate assignment for 14 N gs and the initial mirror T = 1 ground states in the p -2 basis space (even with 
mixing) has no L — 2 component. Only by mixing in L = 1 and L = states in the T = sector can the mirror ground 
states have non-zero MGT elements. In the shell model the use of reasonably strong tensor and spin-orbit interactions 
causes the L = 2 state in 14 N gs to become significantly mixed with the L — and L = 1 basis states which has to be 
fine tuned to cancel the MGT elements when the T = 1, L = and L = 1 states in 14 C and 14 are automatically 
strongly mixed by the strong spin-orbit interaction. The nature of the reference 16 state in this alternative model 
leads to the expectation of a strong suppression of the nuclear spin-dependent mixing interactions which is why the 
model has weak mixing between the LSJ states for both T = and T = 1 states and changes the infinite life of 14 C 
into a long lifetime. The model explains the large asymmetry between the mirror /3-decays because of the interference 
between the nuclear and electromagnetic spin-orbit mixing interactions with the latter term having opposite sign for 
the state in 14 to that in 14 C. The MGT elements for these two transitions are indeed very different in magnitude 
with this model and also their MGT elements have opposite sign in order to obtain the observed negative values for 
the a+, a_ shape slope parameters. Hopefully the detailed results provided in this work on the /3-decays in A = 14 
nuclei will be helpful to future investigators in their search for a more complete description of the structure of p-shell 
nuclei. 



Acknowledgments 

My colleagues at Florida State University in particular K. Kemper, A. Volya, V. Abramkina, and S. Tabor are 
greatly acknowledged. Special thanks for several discussions are due to J. C. Hardy, B. A. Brown, A. Richter and D. 
J. Millener. 



[1] D.R. Inglis, Rev. Mod. Phys. 25, 390 (1953). 

[2] B. Jancovich and I Talmi, Phys. Rev. 95, 289 (1954). 

[3] W.M. Visscher and R.A. Ferrell, Phys.Rev. 107, 781 (1957). 

[4] L. Zamick, D.C. Zheng and M.O. Fayache, Phys. Rev. C51, 1253 (1995). 

[5] M.S. Fayache, L. Zamick and H. Muther, Phys. Rev. C60, 067305 (1999) 

[6] R. Sherr et al, Phys.Rev. 100, 945 (1955). 

[7] J. P. Elliott, Phil. Mag. 1, 503 (1956). 



14 



[8] I. Talmi, Phys. Rev. C39, 284 (1989). 

[9] S. Aroua et al, Nucl.Phys. A720, 71 (2003). 
[10] J.W. Holt et al, Phys. Rev. Lett. 100, 062501 (2008), (see also J.W. Holt, N. Kaiser, and W. Weise, Phys. Rev. C79, 

054331 (2009)). 
[11] H.Genz et al, Z. Phys. A341, 9 (1991). 
[12] R.L. Huffman et al, Phys.Rev. C35, 1 (1987). 

[13] K. Amos, D. Koetsier and D. Kurath, Phys. Rev. C40, 374 (1989). 
[14] J.A. Wheeler, Phys. Rev. 52, 1083 (1937). 

[15] D.M. Dennison, Phys. Rev. 57, 454 (1940), ibid 96, 374 (1954). 

[16] S.L. Kameny, Phys. Rev. 103, 358 (1956). 

[17] M. Freer, Rep. Prog. Phys. 70, 2149 (2007). 

[18] D. Robson, Nucl. Phys. A308, 381 (1978). 

[19] D. Robson, Prog, in Part, and Nucl. Phys. 8, 257 (1982). 

[20] D. Robson, Phys. Rev. D35, 1018,1029 (1987) and Phys. Rev. C61, 015202 (1999). 
[21] C. Long and D. Robson, Nucl. Phys. A474, 627 (1987). 

[22] J.S. M c Carthy, I. Sick and R.R. Whitney, Phys. Rev. C15, 1396 (1977), (see also R. Arnold et al, Phys. Rev. Lett. 40, 
1029 (1978)). 

[23] J.A. Nolen and J. P. Schiffer, Ann. Rev. Nucl. Sci. 19, 471 (1969). 

[24] I.S. Towner and J.C. Hardy, Phys. Rev. C72, 055501 (2005) and J. C. Hardy, private communication (2010). 
[25] W.T. Chou, E.K. Warburton and B.A. Brown, Phys. Rev. C47, 163 (1993). 
[26] A. Garcia and B.A. Brown, Phys. Rev. C52, 3416 (1995). 
[27] F.E. Wietfeldt et al, Phys. Rev. C52, 1028 (1995). 

[28] G.S. Sidhu and J.B. Gerhart, Phys. Rev. 148, 1024 (1966), (see also ref. [H for corrections to C(Z,W)). 

[29] E.K. Warburton and B.A. Brown, Phys. Rev. C46, 923 (1992). 

[30] A. Negret et al, Phys. Rev. Lett. 97, 062502 (2006). 

[31] H. Crannell et al, Nucl. Phys. A278, 253 (1977). 

[32] R.S. Willey, Nucl. Phys. 40, 529 (1963). 

[33] R.P. Singhal et al, Phys. Rev. C28, 513 (1983). 

[34] R.L. Huffman et al, Phys. Lett. 139B, 249 (1984). 

[35] N.J. Stone, At. Data Nucl. Data Tables 90, 75 (2005). 

[36] J.W. de Vries et al, Phys. Lett. 205B, 22 (1988). 

[37] F. Ajzenburg-Selove, Nucl. Phys. A523, Table 2 (1991). 



